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ABSTRACT 



We measure the broad emission line region (BLR) size of a luminous, L ~ 10 ergs , high-z quasar using 
broadband photometric reverberation mapping. To this end, we analyze ^7.5 years of photometric data for 
MACHO 13.6805.324 (z ~ 1.72) in the B and R MACHO bands and find a time delay of 180 ±40 days in 
the rest frame of the object. Given the spectral-variability properties of high-z quasars, we associate this lag 
with the rest-UV iron emission blends. Our findings are consistent with a simple extrapolation of the BLR 
size-luminosity relation in local active galactic nuclei to the more luminous, high-z quasar population. Long- 
term spectroscopic monitoring of MACHO 13.6805.324 may be able to directly measure the line-to-continuum 
time-delay and test our findings. 

Subject headings: galaxies: active — methods: data analysis — quasars: emission lines — quasars: individual 
(MACHO 13.6805.324)— techniques: photometric 
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1. INTRODUCTION 

Our understanding of galaxy and black hole (BH) forma- 
tion a nd co-evolution has progressed significantly in recent 
years (IWandellll999tlFerrarese & Merrittll2000| lHaring & Rixl 
2004t IPeng et alj 120061: iBennert et all 120101: Icisternas et all 



201 II) with new observations being able to probe the first 



epoch of quasar activity, and place interesting constraints on 
the me chanisms responsible for BH growth and galaxy for - 
mation dHughes & Blandfordl200'3l:lTrakhtenbrot et alJ201 lb . 
Gaining further physical insight requires that the BH masses 
be determined with good accuracy in a large sample of ob- 
jects. 

Presently, the best means for weighing BHs in quasars is vi a 
the reverberation mapping (RM) technique (lPetersonlll993l) . 
which measures the size of the broad line region (BLR). This, 
combined with a measure of the velocity dis persion of the 
BLR, can be used to estimate the BH mass ( Peters on et alj 
2004). To date, BH masses in ~ 45 low-z objects have been 
measured in this way leading to various scaling laws with 
other quasar properties. These relations are often extrapo lated 
to hig h-z objects to allow for the indirect estimate of their BH 
mass (Netzer 2003). Nevertheless, it is not clear that such ex- 
trapolations are meaningful, and a more direct measure of the 
BH mass in such objects is highly desirable. 

Here we report new measurements for the BLR size 
in a z — 1.72 luminous quasar, using photometric RM 
dChelouche & Da niel 2012|), an d by analyzing ^ 7.5 years of 
data from the MACHO survey dGeha et alj|2003l) . Section 2 
presents the analysis and results, with a follow-up discussion 
in §3. 

2. ANALYSIS & RESULTS 

To date, some 200 quasars have been confirmed behind the 
Magellanic clouds dDobrzvcki et al. l l2002t iGeha et all 120031: 
lKozIowskietaU2012l) . In the course of analyzing their photo- 
metric data, we report our findings for MACHO 13.6805.324, 
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an niR = 18.66 quasar at z — 1.72 having a monochromatic 
luminosity AL A (1350A) - 3 x ltf^ergs -1 . We focus on 
this object as it has the best-sampled B a nd R light curves, 
as well as the highest fractional variability ( Kasp i et~aLll2007L 
and references the rein), F yw ~ 0.16, of a high-z quasar in the 
IGeha et all d2003l) sample. This object is consistent with be- 
ing radio-quiet down to a flux limit of ~ 6 mJy ( 5 mJy) at 
21 cm (13 cm) dMarx et alll997tlMauch et al.ll2003l) . A more 
complete analysis of the MACHO quasar sample is deferred 
to a forthcoming publication. 

We choose to work with the original MACHO bandpasses 
rather than transform into standard B and R bands. This 
reduces the spectral overlap as the photometric transforma- 
tions require linear combinations of the MACHO magnitudes. 
Light curves were filtered against poor data: all points with er- 
rors larger than £ times the mean error or those which deviate 
from the light curve mean by more than £ standard deviations 
were discarded (Fig. [T}. 

2.1. Spectral decomposition 

To effectively use the photometric RM technique of 
IChelouche & Daniel d2012l) . it is helpful to identify the band 
with the larger contribution of emission lines to its variance. 
If the emission line contribution to the flux may be used as 
a proxy to its variable component (e.g., as would be the case 
if the relative flux variations in all emission lines were sim- 
ilar; see §3), then by spectral decomposition and knowledge 
of the instrumental throughput, one can identify the line-rich 
and line-poor bands. Unfortunately, the quality of available 
spectra of MACHO 13.6805.324 is low, and the data are not 

4 To estimate the monochromatic luminosity of MACHO 13.6805.324 
at 1350A, we take the reported R-magnitudes from IGeha et all 
12003), employ a standard K-correction, an d account for the e xtinc- 
tion behind the LMC using the results of Zaritskv et al] <2004l see: 
http://ngala.as.arizona.edu/dennis/lmcext.html I. To this end, we average over 
50 stars of all types in the direction of the quasar, and identify the upper 
turn-over in the extinction distribution, at Ay ^ 1.3 mag, with the effective 
extinction of the background source. We conservatively estimate the uncer- 
tainty on the luminosity to be a factor 3, to account for the potentially patchy 
structure of LMC disk. Our estimate is consistent with mg-mn for MA- 
CHO 13.6805.324 being larger by ~ 0.67 mag than typical, as concluded by 
integrating the Vanden Berk et al. (2001) composite quasar spectrum over the 
MACHO bands, and assuming a Galactic extinction law. Using concordance 
(0.3,0.7,0.3) cosmology, we find AL^(1350A) ~ 3 x 10 46 ergs-'. 
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Figure 1. R (upper panel) and B (lower panel) light curves for MA- 
CHO 13.6805.324 cover ~ 7.5 years in the observed frame (~ 2.8 years in 
the quasar frame). The number of points in each band is of order 10 3 , and 
depends on the filtering scheme used: darker shades correspond to effectively 
unfiltered data (£ = 10) while brighter colors use a f = 2 filter. The fractional 
variability measure, F var ~ 0.14 (0.16) when £ = 2 (£ = 10) filter is used. 
The mean photometric error is ~ 5.5% (7.6)% for £, = 2 = 10) . 

flux-calibrated. Nevertheless, t heir broad consistency with the 
composite quasar spectrum of iVanden Berk et"aTI (120011 see 
Fig. [2]) motivates us to use the latter for our purpose. 

Prominent emission lines at z = 1.72 include the 
C III] A 1909, Mg IIA2799, as well as the iron complexes 
(Fig. |2). To account for the iron-blen d contribution to the 
spectru m, we use the iron template of Vestergaar d & Wilkesl 
(2001). An underlying powerlaw continuum model was as- 
sumed, Fx A.~ 1,63 , which is somewh at different from the 
one used by IVanden Berk et al.l (|2001), and better accounts 
for the flux level at ~ 6000A (observed), where little iron 
blend and Balmer continuum emission is present. Emission 
lines and blends were individually convolved with a single 
Gaussian kernel of some amplitud e and width, and a qualita- 
tive agreement was sought with the Vanden Berk et al. (2001) 
composite spectrum (Fig. [2j. The solution is not strictly 
unique, nor very physical, and its sole purpose is to qualita- 
tively assess the relative contribution of the various emission 
components to the broadband flux. 

Weighing the decomposed s pectrum by the instrumental 
throughput dAlcock et al.ll 1999b . we find that the fi-band has 
the greater relative contribution of emission lines to its flux. 
Specifically, the iron emission blends contribute a net of ~ 
10% to the flux in the bands, with the relative contribution 
of the Mg II 12799 and C III] A 1909 to their respective bands 
being > 4 times smaller. 

2.2. Uncovering the time-lag 

Having identified the line-rich band, we now proceed 
to measure the line - to-c ontinuum time-delay. Following 
IChelouche & Da niel (2012|), we compute two statistical es- 
timators for the line-to-continuum cross-correlation function: 
€ca(5*) = /r * h ~fB* /b, and €aa(5*) = f R */« - f B * fy 
(" * " denotes convolution and fs, /r are the light curves in 
the B and R bands, respectively; see Fig. |3). With these defi- 
nitions, a peak at 8t > indicates an emission component in 
the B-band, which lags behind the B-band. To calculate these 
statistical estimators, we use two indep endent schemes: the 
interpolated cross-correlation function (Peterson et al.ll2004l 
ICCF) and the z - transformed discrete correlation function 
( I Alexander! 1 1 9971 ZDCF). The former has the advantage of 
being somewhat more sensitive, while the latter is less af- 
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Figure 2. A quasar composite spectrum at z = 1.72 fVanden Berk et all 
120011) overlaid with the throughput curves of the MACHO B- and R-filters 
[the somewhat erratic behavior is real and results from an interference pat- 
tern (C. Stubbs; private communication)]. A qualitative spectral decompo- 
sition into the prominent emission components is shown including the iron 
blend (gray curve), the C III] A 1909, and Mg II A 2799 lines (magenta curves), 
and a powerlaw continuum (dashed line). Designated minor emission lines 
were also included in the spectral decomposition (green curve). Integrat- 
ing over the emission line flux, and taking into account the wavelength- 
dependent throughput, the largest contribution to line emission is from the 
iron complexes in the R— band. The contribution of the C III and the Mg II 
lines is roughly 4 times smaller than iron's and their relative contribution 
to their respective bands is similar. Also shown are the two {rum flux- 
calibrated, slightly shifted in wavelength range for clarity) optical spectra 
for MACHO 13.6805.324, which were multiplied by a power-law function of 
wavelength so that their trend roughly follows that of the composite spectrum 
(dashed blue-shaded curves). The spectral shape is qualitatively consistent 
with the quasar template (note the instrumental feature at around 7300A) and 
suggests a significant contribution of the iron emission blend to the flux in the 
S-band. Note the prominent atmospheric absorption feature which happens 
to coincide with the Mg II emission line at this particular redshift. 

fected by sampling. 

Results for <^ca an d <Haa using the ICCF and ZDCF 
schemes, and for two levels of light curve filtering, are pre- 
sented in figure [3] All analyses give consistent results: (*ca 
and <^aa peak at 8t ~ 400 — 600 days. Unsurprisingly, the 
peak is more significant when modest filtering (£ = 2) is 
applied, resulting in £,ca, <^aa being significant at the ~ 
2.8(7, 5(7 ((7 is the standard deviation) level, respectively. 

To verify that the observed signal is indeed associated with 
a lagging component in the R -band, we reversed our choice 
of line-rich and line-poor bands, and repeated the analysis: 
no significant peak was detected at 8t > (instead, a highly 
significant trough was detected at St ~ 400 — 600 days; not 
shown). In addition, we carried out sets of Monte Carlo 
simulations (Chelouche & Daniell l2012) wherein the Z?-band 
light curve was convolved with a Gaussian kernel having pre- 
scribed time-delay and width (the latter was taken to be half of 
the former but the results are not very sensitive to the particu- 
lar kernel chosen). The resulting light curve was then scaled 
down and added to the B-band signal so that the line contribu- 
tion to the flux is ~ 10%. The combined light curve was then 
sampled with the cadence of the B-band to create a synthetic 
light curve, The ^-estimators were calculated for the light 
curve pair [/b,/^], and the results for several input time-lags 
are shown in figure [3] Clearly, the observed signal is qualita- 
tively reproduced by our simulations for lags of the order of 
400 days. This a) confirms that the observed signal is indeed 
associated with a lagging emission component in the /?-band 
and is not spurious, and b) that the time-lag constraints de- 
duced below are meaningful and are not significantly biased. 

We estimate the time-lag and its uncertainty in the fol- 
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Figure 3. Photometric RM analysis. Top panel: an excess power is seen at around 500 days in the cross-correlation signal compared to the auto-correlation 
signal (thick lines). The simulated signal qualitatively reproduces the excess power in the cross-correlation term, at the right timescales, for three input time lags 
(see text and legend). Bottom panels: the calculated ^-functions for two levels of filtering: £ = 2 (left) and £ = 10 (right). Solid curves show the results of the 
ICCF algorithm with shaded regions indicating the uncertainty. The ZDCF calculations are shown as points with error bars in dotted lines. Blue (red) colors 
correspond to ijcA (£aa)- All methods give consistent results and a modest level of light curve filtering yields a more significant result. We constrain the time 
scale at which the £ -functions peak to be ~ 500 ± 100 days in the observed frame. Repeating the analysis using model-dependent synthetic R-band light curves 
(blue-shaded curves; see text) shows an agreement with a time lag of >400 days; considerably shorter or longer lags are excluded by the data. 

lowing way: the peak in 4ca an d 4aa is identified with the 
time-lag for each Monte Carlo realization (a total of 100 re- 
alizations are used for each numerical scheme). A time- 
lag distribution is then obtained whose mean is identified 
with the time lag and its uncertainty with the standard de- 
viation. We find that the time-lag is 480 ± 100/530 ±70 days 
(ZDCF), or 490 ± 80/490 ± 50days (ICCF) for the ^ C a/4aa 
statistical estimator. We conservatively estimate the lag to 
be 500 ± 100 days, which corresponds to a rest-frame lag, 
X~ 180±40days. 

3. DISCUSSION & CONCLUSIONS 

With only photometric data at hand, interpreting our find- 
ings is not straightforward: while we have a statistically ro- 
bust time lag measurement, it is difficult to associate it with 
a specific emission linfl That being said, it is very likely 
that the observed time-lag is associated with the iron emis- 
sion blends for the following reasons: 1) a positive peak in 
4cA ; 4aa indicates that emission lines contributing to the R- 
band are responsible for the signal, 2) the iron contribution to 
the /?-band typically exceeds that of the Mg II 2799A line by a 
factor > 4 (Fig. |2), and 3) the implied flux variatio n of the iron 
blend is consistent with that seen in l ow-z objects (Maoz et al. 
Il993t IVestergaard & Peterson 2005), and indirectly deduced 
for high-z quasars (Meusinger et al. 201 10 

5 Another possibility is that the signal originates in an unidentified contin- 
uum emission component. While this cannot be excluded based on available 
data, we find it rather unlikely given that continuum emission, by definition, 
extends over a broad wavelength range, and the fact that there is no clear 
reason why the lag should be similar to that which characterizes emission 
lines. 

6 In a recent work, Meusinger et al. (2011) analyzed Stripe 82 SDSS 
quasars and found that the combined flux from all emission lines and blends 
varies by ~ 10%, which is consistent with the relative variations measured 
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Figure 4. The arbi trarily normalize d (non flux-calibrated) spectra of MA- 
CHO 13.6805.324 iGeha et al. 2003, and S. Kozlowski, private communica- 
tion; see legend) show two prominent emission lines. A cubic spline was fit 
far from the prominent line locations to normalize the spectrum (not shown), 
and single Gaussian models were overlaid to constrain the FWHM range, as 
indicated by the shaded areas. Only the red wing of the C IV A 1549 line was 
used for fitting purposes due to potential absorption features just blue-ward 
of its peak (left panel). Similarly, only the red wing of the C III] A 1909 line 
was used for FWHM estimation to avoid blending with other emission lines 
just blueward of its center (right panel). We note that FWHM measurements 
for this line should be treated with caution since iron blend emission may be 
substantial (see also Fig. ff), and is dififcult to estimate given the poor S/N of 
the non flux-calibrated data. 

Little is known with confidence about the physical prop- 
erties of the UV iron-emitting region, and a statistically sig- 
nificant time-lag, reflecting on its size, exi sts for only one ad - 
ditional low-luminosity source, NGC5548 dMaoz etal.11 993). 



for individual lines by Kaspi et al. 1 2007). This implies that si milar flux vari- 
ations of the iron blend are characteristics of quasars [see also Wil hite et all 
12005) for some additional constraints on the time variability of emission 
lines in high-z quasars]. 
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Figure 5. The BLR size-luminosity relation (left; adopted from Kaspi et al. 12007); see their paper for a discussion of the va ri ous fitted relations) and the mass- 
luminosity relation [right; adopted from Peterson et al. (2004, and augmented by up-tp-date data from Peterson et al. (2005); Denney et al. (2006); Kaspi et al. 
)2007l) ; lGner et all (2008); Ben tz et all 420091) ; IDennev et alj (2010); Gri eret all ( 120121) 1 using the composite spectrum of Vanden Berk et al. ( 200JJ) to convert 
monochromatic luminosities to AL^ (1350A); uncertainties are not shown]. Left: plotted are data from Peterson et al. (2004, 2005); Kaspi et al. (2005, 2007) as 
well as our res ults for MAC HO 13.6805.324, which are in agreement with both relati ons. Also shown, for comparison, are reliable BLR size measurements for 
the iron blend iMaoz et al. 1993) and the C III A 1909 and Mg II A2799 emission lines I Metzroth et al. 2006); see legend. Evidently, the BLR size measured here 
lies on the extrapolation of the R — L relation, as obtained for low luminosity, low-z objects. Right: Interestingly, MACHO 13.6805.324 is consistent with emitting 
at > 10% of its Eddington rate, in agreement with the efficiency of other bright quasars [note t hat the mass reported here for S5 0836+71 is ~ 1.8 times larger than 
the value reported in Kaspi et al. (20071) an d follows from the normalization of Onken et al. (2004)]. Curves of constant L/LEdd. are also plotted (dashed lines) 
assuming the bolometric correction used by Kaspi et al. (2007). 



Taken together, the results imply that the size-luminosity scal- 
ing for this region is con sistent with that of other broad lines 
(Kaspi et al. 2005L I2007L see our Fig. [5]). Interestingly, the in- 
ferred sizes of the iron- an d C IV A 1549-em itting regions are 
comparable in NGC 5548 dMaoz et al.lll993l) . If true also for 
MACHO 13.6805.324 then we find consistency with a simple 
extrapolation of size-luminosity relation for the C IV-emitting 
region, as obtained for lower luminosity sources, to the more 
luminous quasar population (Fig. [5). A meaningful compari- 
son with the size-luminosity relations for other rest-UV lines 
is currently limited by small number statistics. 

Once the size of the line-emitting region is known, the 
black hole mass may be determined via Mbh — 2.7 x 
10 5 (FWHM/10 3 kms- 1 ) 2 (T/da ys)M^ dPeterson et al.l l2004l 
and using the normalization of Onke n et al. I (120041) 1 where 
FWHM is the full width at half maximum of the relevant 
emission line. This expression assumes that the gas is viri- 
alized and uses the FWHM as a proxy for its velocity dis- 
persion in the black hole's potential well. Nevertheless, the 
FWHM of the iron blend in MACHO 13.6805.324 cannot be 
determined from the available spectra, and it is not clear tha t 
the iron-emitting gas is virialized in quasars dHuet al.i r2008). 
With these uncertainties in mind, we note that currently avail- 
able data point to a roughly comparable size f or the iron, 
C IV A 1549, and C III A 19 09 emitting regions (IMaoz et alj 
[1991 lMetzrothetaLll2006l see our Fig. EE If true also 
for MACHO 13.6805.324 then we may use the FWHM of the 
carbon lines and our iron line region size measurement to es- 
timate MbiE Taking FWHM=6500± 250010ns- 1 (Fig. g), 



we estimate the black hole mass to be 2 x 10 9 M© with a for- 
mal uncertainty of ~ 0.3 dex (see Fig. q). We note, however, 
that this mass estimate does not account for potentially impor- 
tant, but poorly understood, systematic effects concerning the 
physics and geometry of the iron-blend emitting gas. 

It is worth noti ng that, being > 10 quieter in the radio band 
than S5 0836+71 kaspi et alJ2007h . MACHO 13.6805.324 is 
not a blazar, and may well be a radio-quiet object. As such, it 
may have the largest black hole mass ever measured using the 
RM technique in a representative member of the quasar pop- 
ulation. Our findings also indicate that MACHO 13.6805.324 
is probably shining at > 10% of its Eddington luminosity 
(Fig. |5]>, and in agre ement with similar quasars in its class 
jShemmer et al.ll2006l) . 

A more reliable interpretation of our results, as well 
as a more accurate BH mass determination, require a 
higher signal-to-noise spectrum. A sparse (one visit ev- 
ery two months) spectroscopic-monitoring campaign of MA- 
CHO 13.6805.324, over a period of - 7 years, will deter- 
mine the variability amplitude of the iron blend as well as po- 
tentially (independently) measure the time-lag in this object. 
Furthermore, by analyzing large enough samples of MACHO 
quasars, with a similar redshift and luminosity range to MA- 
CHO 13.6805.324, it may be possi ble to statistically corrob - 
orate the results presented here (IChelouche & Daniel] 120 121) . 
With these future tests and improvements in mind, our present 
findings already demonstrate the feasibility of long-term pho- 
tometric surveys in determining the BLR size (and BH mass) 
in luminous high-z quasars. 



The situation concerning the size of the C III A 1909 region is far from 
being clear: two low-luminosity objects having only loose constraints on 
the ti me-lag indicate that it may be twice the size of the C IV A 1 549 re- 
gion {Peterson 1994j), while a tentative time-lag measurement in a high-z 
quasa r indicates that th e line traces the continuum level with no apparent 
delay IKaspi et all2W7l) . 

8 There are uncertaintie s associated with the use of carbon lines for mea- 
suring M B h : as noted by Baskin & Laor ( 2()05|), the C IV A 1 549 FWHM may 
be affected by non-virial motions, especially for bright objects emitting close 
to their Eddington limit. Also, line de-blending may be required to prop- 
erly estimate the FWHM of the C III A 1909 line. Lastly, a generic problem 
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with single-epoch FWHM estimates concerns the fact that the FWHM of the 
varying comp onent of the line may differ from that measured from the mean 
spectrum ( Kaspi et al. 2000). 
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